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PREFACE 



This paper presents the results of a study of Constant 
Absolute Vorticity Trajectories applied to isentropic surfaces. 
The objectives of this study were; first, to show that Constant 
Absolute Vorticity Trajectories are theoretically better applied 
to forecasting future positions of parcels on isentropic surfaces 
than on constant pressure surfaces as they are now applied; sec- 
ondly, to show quantitatively the actual improven:ent of the tra- 
jectories on isentropic surfaces by conparing statistically 
results of actual forecasts from both types of charts; thirdly, 
to suggest how the improved technique can aid forecasting, par- 
ticularily upper air prognosis. 

This work was undertaken as the thesis requirement for the 
degree of Vaster of Science in Aerology at the U. S. Eaval Post- 
graduate School, Vonterey, California, during the academic 
year 1952-1953. 

The author is indebted to Associate Professor George J. 
Galtiner of the Aerology Departiient for his very valuable assis- 
tance and constructive criticism diuring the investigation. l!e 
also wishes to aclmowledge the assistance of Professor A. Boyd 
Ifewborn of the ^'athe^atics Department and the author’s \d.fe, 
Pauline, who contributed so much help with the laborious task 
of decoding, computing, and entering isentropic charts. 
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I. If’TilOaUCTIOiJ 



The meteorologist is often called upon to forecast for periods 
in excess of tiventy-four hours for which sir-pie extrapolation is 
generally insufficient. One method of forecasting entails prepa- 
ration of prognostic surface charts and preparation of forecasts 
from these charts. As most meteorlogists loiow, making accurate 
prognostic surface charts is easier said than done. One of the 
important tools used in estimating the position and intensity of 
surface systems is a prognostic 500 millibar chart. At present ^ 
much attention is being paid to the problem of upper air prognosis, 
and considerable information has been published concerning fore- 
casting the 500 millibar sxorface, the most recent being ^orecastijig 
in the Middle Latitudes by Riehl and collaborators . It is hoped 
that this paper may shed some light on a facet of this problem and 
nay contribute to better prognosis and forecasting. 

One aid to prognosis of upper air charts is the construction 
of forecast air parcel trajectories based on the principle of con- 
servatism of the vertical component of absolute vorticity,^^ +<f ^ , 
where ^ is the vertical component of relative vorticity and *f 
is the vertical component of the coriolis force loiowi as the 
Coriolis parameter, lossby W initiated the vorticity concept 
by showing that under certain restrictive assumptions, the follow- 
ing relationship holds; 
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The assumptions referred to above are as follows: 

a. The atmosphere is barotropic. 

b. The atmosphere is a homogeneous, incompressible fluid. 

c. Motion is purely horizontal. 

d. Friction forces are neglected. 

e. There is no horizontail divergence. 

Rossby ^9, pp 268-289J then integrated his fundamental equation 
and was able to obtain trajectories for particles which conserve 
their vertical component of absolute vorticity. This trajectory 
is Icnown as the Constant Vertical Component of Absolute Vorticity 
Trajectory and \d.ll be referred to as a CAV trajectory. 

The assumptions necessary to derive the vorticity equation 
according to Rossby are obviously quite restrictive. Starr, as 
editor of the Journal of Meteorology (Jvuie, 1945), suggested that 
CAV trajectories might better be depicted on an isentropic chart 
than on a constant level chart as there are no solenoids on an 
isentropic chart. Namias js, pp 372-374J states that from obser- 
vations available isentropic surfaces appear, to the first approxi- 
mation, to be substantial surfaces, i.e., svu'faces that contain the 
same air particles from day to day. Influences of non-adiabatic 
nature are appreciable over long periods of time, but are usually 
insufficient to disrupt the fundamental isentropic flow from one 
day to the next. The derivation of the vorticity equation for flow 
on an isentropic surface requires fewer assvimptions . Such restrictions 
as purely horizontal flow, a barotropic attnosphere, and a homogeneous 
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incompressible atmosphere nay be renoved. Now, an isentropic 
surface is a surface of constant potential temperature. In 
meteorlogical processes, potential temperature is one of the more 
conservative elements by which parcels of air nay be identified 
from time to time. Since the atmosphere is normally stable, 
potential temperature increases \vLth height cind the atmosphere 
may be considered to consist of an infinite number of isentropic 
surfaces. Sir Napier Shaw JlO, p 26^ first suggested that charts 
of isentropic siurfaces should be dra^vn as the motion of air is best 
resolved on isentropic surfaces. It i^ill be shown in this paper 
that the prognostic trajectory based on the conservatism of the 
vertical component of absolute vorticity is theoretically and 
practically better suited to use on isentropic charts than on 
constant pressure charts as it is now applied. 
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II. THEORETIC.\L II^TISTIGATION 



The motion of particles on a substantial surface can be 
described with Lagrangian methods. The form of the vorticity 
equation for flow on an isentropic surface derived here using the 
Lagrangian method of attack is due to Raltiner''*’. LetX,^^Z and-fc 
be normal Cartesian coordinates and let]X)Y^ ^ r be 

Lagrangian coordinates where 'X and X 3,re coordinates of the 
horizontal projection of the isentropic surface, i. ~ ^ and 'Y 
refers to time for the isentropic sitrface. Then 

chc 






= ^ 
C>^ 



— Q.P c) 









It follows from the above equations that 
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where r=c^T + «^«. 



the stream function for an isentropic surface. 



In a sir.dlar way we find 
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Associate Professor C. J. Kaltiner, U. S. Kaval Postgraduate School 
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Further, since the motion is assumed to be isentropic, 
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Differentiating (7) with respect to X and (8) ivith respect to 
and subtracting, we obtain 

^ I -f. Ou. g)u. a ^ <pt>* <Jfa- 1 XT ' c)^«C _ 
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and introducing vector notation ve obtain 



- (5^ W K.) ~x ^ ^ 

iidiere the operator ^ ^ ^ 1 n ^ 

a>X '^<3X ’ 

Finally, we may wite this as 

since'V — t for any particle. NoX'/, assurae the divergence tem^ .V=0. 
Then we may wite the vorticity equation in the familiar forr.i 




V/hile this form is similar to the form derived by Rossby, there 
are some fundamental differences. Parcel notion is not limited 
to the horizontal, but is tliree dimensional. Moreover, a barotropic 
and incompressible atmosphere is not required. 
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III. TECIKIQUE OF INVESTIGATION 



The relative vorticity may be written in a fom similar to 
the usual expression for the vertical component of vorticity in 
natvural coordinates 

y ^ V ^ 

where P and f\J are measured on the isentropic map. If we select 
a parcel in a broad uniform flow or in the axis of a jet stream, =0. 

Further, if the parcel is at an inflection point in the flow, rr O. 
Under these conditions the CAV trajectory can be computed. It 
will be noted that defined above is not exactly equal to the 
projected radius of curvature of the streamline on the isentropic 
sxirface. However, the approximation involved here is normally 
less than those made in the actual computation of CAV trajectories 
as described by Fultz [4, p 13] and revised by Wobus [12] . hliile 
approximations have to be made, it must be remembered that the true 
worth of a technique lies in the statistically proven results 
obtained. 

The simplest method of obtaining a CAV trajectory is through 
the use of a mechanical device such as the so-called 'Nv'iggle-V7agon'* 
devised by Wobus [|l2j and cvurrently in use in the h’BAN Analysis 
Center, V/ashington, D. C. Original studies of the CAV trajectory 
assumed a plane earth. This assumption entailed errors due to 
distortion when compared to trajectories along a spherical earth. 

The Wobus '^Wiggle-Wagon” corrects this defect. To aid meteorolo- 
gists who do not have access to the "’Viggle-Wagon”, the U. S. Navy 
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Bureau of Aeronautics Project ARO^wV Lu] has published a table 
for conputing CAV trajectories r,iadc from trajectories traced by 
the **Wiggle-yagon" . This table ^vlll be used to compute CAV tra- 
jectories in this paper. For application, it is necessary to Icnow 
speed, direction, and the latitude of the parcel, all at an inflection 
point in the flow. 

The next step in the investigation was to compute trajectories 
on isentropic and 500 millibar charts and test the results statisti- 
cally. A rather complete compilation of upper air and surface data 
are available in the Northern Hei;iisphere Surface and 500 Nillibar 
Charts Series published by the U. S. V/eather Bureau. These remarkable 
publications not only have a daily series of logically analysed surface 
and 500 millibar charts, but also contain a compilation of all data 
used to plot these charts. Data were taken from the period of 
January-February 1949. From this data, a set of thirty isentropic 
charts was plotted. The 303K isentropic surface was selected because; 
(1) this surface was often close to the so-called "level of non- 
divergence" described by Bjerknes [1] ; (2) this surface often con- 
tained much of the 500 millibar ^>rind speed maximum axis (jet); 

(3) the 303K isentropic surface was usually sufficiently far from 
the surface to remove effects of siurface friction and turbulence. 

It was desired to select an isentropic surface containing much of 
the 500 inillibar isotach maxirAun axis so computation of CAV tra- 
jectories could aid in the prognosis of the 500 millibar surface. 
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Inasnuch as the derivation of CAV trajectories on an isentropic 
svirfacc assumes the divergence tcn:i to be aero, the selection of 
an isentropic surface close to the level of no horizontal divergence 
r.iay be desirable. The effect of the divergence tern is discussed in 
Chapter IV. 

After the heights of the 303II isentropic surface were plotted, 
all available wind data pertinent to the surface v;ere plotted. Vdiere 
^vind observations v/ere available from stations other than radiosonde 
stations, an interpolated value of the height of the isentropic surface 
and the wind nearest that height were entered. Because of the fairly 
dense nettvorlc of upper air word observing stations in the United 
States, sufficiently accxirate streamline analysis was possible \d.th- 
out the tir:e-consuning computations of stream function values. 

Streeimline analyses were made for 0300Z each day in the period 
from 6 January, 1949, to 22 Januan'’> 1949, and from 10 February, 1949, 
to 20 February, 1949. For each chart, a nvtnber of CAV trajectories 
were computed. The following conditions had to be met before choosing 
the initial point; (1) the point had to be in or very near an inflection 
point in the flow (undergoing no effects of curvature); (2) the point 
had to be in or near the axis of the isotach maxiniun (jet stream); 

(3) the point must be at an actual point of id-nd observation to obtain 
correct initial wind speed and direction. This procedure netted from 
one to foiu' points suitable for computation of CAV trajectories per 
chart. 
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For comparison with each CAV trajectory plotted on the 
isentropic chart, a companion CAV trajectory was computed on the 
500 millibar chart of the sane time. It was not too difficult 
to find a companion inflection point on the 500 millibar chart as 
the isentropic surface was usually close to the 500 millibar surface. 
After the CAV trajectories were constructed on each chart, forecast 
;djnd velocities were determined by moving the parcel along the tra- 
jectory a distance equal to the initial vrind speed multiplied by 
the number of hours to forecast time. The forecast times selected 
were twenty-four, forty-eight, and seventy-two hours. These fore- 
cast winds were then plotted on the verifying charts and verified. 

On the isentropic chart, verification depended on the forecast 
position falling within the observational network. If the forecast 
position fell within the neti^^ork, but not on a particular station, 
the forecast v^inds were verified by linear interpolation between 
nearby stations. Kov/ever, care was taken to use only stations 
vdiere linear interpolation would be reasonable. On the 500 millibar 
chart, forecast ^vLnds were verified in a similar manner if the fore- 
cast position fell within the observational network or on a reporting 
station. If the forecast position on the 500 millibar chart fell on 
neither of these places, then contours were assicned to be streamlines, 
and gradient ^\dnds were neasiu’ed to verify. Since verification on 
isentropic charts depended on forecast positions falling in the North 
American netA/ork or on a particular station, the number of possible 
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verifications fell off drastically after forty-eis^it hours. However, 
with heiiiispherical 500 millibai' charts available, forecasts rmde on 
the 500 millibar chart could always be verified. All 500 millibar 
CAV trajectories were carried out to seventy-two hours. 

In verifying forecast winds, observed wind direction error 
was arbitrarily selected negative if the observed wind direction was 
to the right (cloclnvise) of the forecast idLnd direction and positive 
if the observed wind direction was to the left (counter-cloc!<wise) 
of the forecast ;d.nd direction. Similarily, wind speed error was 
selected negative if observed wind speed was less than forecast and 
positive if the ^dnd speed was greater than forecast. Ml forecast 
winds and verifications were tabtilated in a fora suitable for 
statistical testing. 

Histograms of wind speed error and wind direction error were 
prepared for twenty-four and forty-eight hour forecast times for both 
the isentropic and 500 millibar charts. These histograms are depicted 
as Figures 1 through 4. Wind speed errors were grouped into cells 
each having a cell interval of 10 Imots error. The cell midpoints 
are for± O, etc., imots error. Since the cell boxmdaries 

are exactly halfway between the cell midpoints, some of the obser- 
vations fell on the cell boundaries. In this case, it was deemed 
best to divide such observations bet^^een each adjacent cell. Less 
trouble was encountered in grouping wind direction error. Ml v\dnd 
direction errors were in incrraents of 10 degrees and values ofO,±!0 
± 20 , etc. degrees of error were taken as cell midpoints. 
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After the data were grouped in histograris, normal curves were 
fitted in accordance with the method described by Hoel ^5, pp 191-194J 
After obtaining these fitted nomal curves,^ tests were performed to 



determine if each sar^iple of data could reasonably be from a nomal 
population. In all cases the sample distributions could reasonably 
be said to be from norraal populations, although in two cases the 



test were met by grouping cells with too few frequencies together. 

After showing normality of the data, the distributions of data from 
the isentropic charts were compared to similar distributions from 
the 500 millibar charts by the use of the test. The test 
v^as made to determine if compared data distributions from the 
isentropic and 500 millibar charts could statistically be from the 
same parent normal population. The purpose here is to show that a 
given distribution of data from the isentropic charts is not from the 
same nomal population as the similar data distribution from the 500 
millibar charts. If there is a significant difference, then inspection 
and intuition v;ill reveal which of the two types of charts are better 
suited to forecasting by means of CAV trajectories. 




value . 
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IV. RESULTS 



It is well to look first at the data in ’’everyday” terns to 
find what kind of results a forecaster nay expect. If it is arbi- 
trarily decided that a wind direction forecast ^<111 verify as a ’’hit” 
if it is \idthin 20 degrees of observed direction, then we find that 
96,4% of the tii?enty-four hour virA direction forecasts from GAV tra- 
jectories on isentropic charts verified as ”hits”. By comparison, 

76.3% of similar forecasts from 500 millibar charts verified as ”hits”. 
V/ith forty-eight hour ^vind direction forecasts, 81, o from isentropic 
charts verified as “hits” and 41% from 500 millibar charts. 

If it is arbitrarily decided that 15 knots of wind speed error 
is not too much to verify as a ”hit”, we find that 76,5 of twenty-four 
hour wind speed forecasts from trajectories on isentropic charts 
verified as "hits”. Of the data from the 500 millibar charts, 42.4,o 
of twenty-four hour wind speed forecasts verified as "hits”. Similarly, 
forty-eight hour wind speed forecasts verified 78.5}5 from isentropic 
chart data and 39,5 from 500 millibar chart data. So, percentage\d.se, 
forecasting \mnd vectors by CAV trajectories on isentropic charts 
was superior to forecasting wind vectors by CAV trajectories on 500 
millibar charts during the period of this study. 

Figures 1 through 4 are histograms of error betw^een observed 
and forecast v:ind speed and direction. Applying standard statistical 
methods, normal curves were fitted to the data. The *y test was 
applied to determine if the distributions of error were likely to be 
from a normal population. Only the twent 3 ’--four and forty-eight houi' 
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errors were compared and tested as the seventy-two hour isentropic 
seunple v/as too small. In each case tested, the distribution was 
shown to be reasonably from a normal population. In cases of data 
from isentropic charts, there was little doubt that the distributions 
were from non'.ial populations. But in cases from the 500 millibar 
charts, the values of*y were close to the critical values and, 
therefore, there is some doubt that they may be from a normal popu- 
lation. However, since the values of y*” did fall below the critical 
values, we will treat such data as being from norr.ial populations. 

The "F" testjs, pp 152-154] is a statistical test for comparing 
variances (>*) of two samples of data to determine if the t^>ro saniples 
can reasonably be from the same normal population. The basic require- 
ment of the "F” test is that the data samples must each be from normal 
populations, although not necessarily the same normal population. This 
is why it was necessary to show normality of the data samples above. 

To apply the "F” test, we first postulate that there is no significant 
difference between samples tested, i.e., = • This is the null 

hypothesis. If the value of F computed falls below a critical value, 
the null hypothesis is not denied and there is no statistical difference 
between the two samples. If the value of F computed is higher than the 
critical value, the null h)rpothesis is denied, and there is a statisti- 
cally significant difference bettv'cen the two samples and they are not 
reasonably from the sam.e parent normal population. It is this last 
condition we wish to obtain here. The results of applying the ”F” 
test are as follows: 
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1. T\/enty-four Hour \7ind Direction Error . The variance of 
data from the isentropic charts 5,’ I.^A . The variance 
of data frori the 500 nillibar charts S^zz^.21 . Then 

r= sT/sT= 

^ at a 95/0 level of belief is 4.90. Therefore, the 
samples tested are significantly different and the null 
hjTJOthesis is denied. 

2. Forty-eight Iloxir >7ind Direction Error . The variance of 

date from the isentropic charts 5,'’-== 3* 1 5 . The 

varicuice of data frori the 500 millibar charts 5^\r 300Z. . 
Then 

F= 5.7 s," = t7S. 

at a 95,0 level of belief is 3.75. This is higlily 
significajit and the null hypothesis is certainly denied. 

3. Ti-renty-four Hour i/ind Speed Error . The variance of data 
from the isentropic charts is 155.7. The variance of data 
from the 500 raillibar charts is 413.1. Then 

F=r 2XS. 

Fc. at a 95,0 level of belief is 5.19 and, therefore, the 
null hj'pothesis is not denied and there is no significant 
difference between the sarnples ef data. 

4. Forty-eight Hour I/ind Speed Error . The variance of data 
from the isentropic charts is 127.3. The variance of data 

from the 500 millibar charts is 515.5. Then 

F= 4.05. 

At a 90, T level of belief the value of F^:r4.00 , and at 
this level the difference bet^/een the saiviples becomes 
significant. 
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Significant results were obtained in the cases of wind direction 
error and in one case of v/ind speed error. In each of the cases 
vdiere significance was obtained, an intuitive ercairdnation of tlie 
data shows that it is the data fx'or. trajectories constructed on 
isentropic charts that are better. 

It is gratifying to note tliat the greatest ii-provencnt was 
showx for xaind direction forecasting for a forty-ed.ght hoxir period. 

In nuch operational forecasting, wind direction is usually deer.’.ed 
nore important than vdnd speed for prognostic work. Vhe proper 
placing of troughs and ridges and associated weather patterns 
depends greatly on correct wind direction forecasts. Therefore, 
from the results obtained, the CAV trajectory method will show 
greatest improvement over the present system of taking CAV tra- 
jectories on constant pressure charts if used on isentropic charts 
to forecast wind direction, at least up to forty-eight hovu’s. 

A recent study of CAV trajectories along a psuedo-600 rdllibar 
chart was published by Bruch Ld . Bruch’s results on wind direction 
error arc quite sii!U.lar to results obtained herein for the 500 rdllibar 
chart. Bruch did not verify xvlnd speeds. 

It is realized that the isentropic chart leaves nuch to be desired 
as an operational weather chart. At present the chart is of course 
laborious to plot since no isentropic data are transratted on tele- 
type circuits; however, if such data were transmitted on teletype 
circuits, strcar.fLine analj’^sis would take little longer to prepare 
than the present 500 millibar chart, .\nother difficulty is that the 
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never. ent of the isentropic surface is z’ot easy to predict, 
the CAV trajectory developed herein is a tv;o dirensional projection 
of a tiu'ee dimensional path, the height of the parcel in the future 
is in doubt. Of course, if the isentropic sxirface did not r.ovc, the 
future position of a parcel would be certain. Out, this is not 
generally the case. A futiu'e study to correlate this study with the 
prognosis of isentropic surfaces could be to advantage. Nonetheless, 
the computation of CAV trajectories on an isentropic surface can be 
of value to prognosis of a nearby constant pressure surface. The 
positions of prognostic troughs and ridges of the isentropic floxv 
can be placed on a nearby constant pressiu'e surface. 

From study of the CAV trajectories on the isentropic charts of 
this paper, it qualitatively appears that parcels ir.’bedded in north- 
westerly flow descend as they move southeastward, whereas parcels 
imbedded in southwesterly flow rise as they move northeasti/ard . 

In the final fom of the vorticity equation there is a diver- 
gence term that is assumed to be zero. This assumption is not always 
fulfilled in the atmosphere. The effect of this divergence temx may 
be important. For illustration, let the divergence term 
be constant in equation (9) and integrating ;d.th respect to -fc , we 
obtain 
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Let*(;=-0 at the initial inflection point, where also Q . 






and 
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As can be seen, positive values of the divergence tern in a north- 
westerly flow \\dll cause a negative contribution to relative 
vorticity of a parcel causing it to go farther sovith than under a 
GAV trajectory. Positive values of the tenn in southwesterly flow 
\ri.ll have a similar effect on the relative vorticity of a parcel 
causing it to not reach the najdLnun latitude indicated by a GAV 
trajectory. Negative values of the divergence term ;/ill cause 
positive contributions to the relative vorticity of a parcel. In 
northwesterly flow this \vlll cause a parcel to curve more cycloni- 
cally than v/ould be indicated by a GAV trajectory and to cui've less 
cyclonically in southv/esterly flow. Fultz |^4, pp 32-10^ mentions 
that in a lower level (10,000 foot level) the average type of 
deviation from GAV trajectories is such that could be caused by 
horizontal velocity divergence in northwesterly flow and convergence 
in southwesterlj'" flow. Mow, 



Thus, the divergence term is composed of t\^ro parts, horizontal 




c>2. ^T/. 



velocity divergence and a shear terxA 



. Fleagle and PanofsI<y 
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have given several values for horizontal velocity divergence, 
with a range of about ^ 5 x 10~^ per second. The values of the 
shear tern could be computed at initial inflection points for 
thirty-four of the CAV trajectories on iscntropic charts of this 
study. In fourteen cases the tem^ was less than 10~^ per second. 

In eleven cases the value of the shear term was approximately 10“^, 
and in nine cases, approximately -10“® per second. Usually, the 
shear term was positive or negative if the flow was up or doivn the 
isentropic surface respectively. From this we see that in about 
60^ of the cases studied, the shear tern was of tlie sane order of 
magnitude as the liorizontal velocity divergence. In these cases 
the shear tern either canceled or doubled the horizontal velocity 
divergence term. Now, if the assir'.ption of V — O were the 

most important assumption in the development of CAV trajectories, we 
night expect forecasting of future positions of parcels by CAV tra- 
jectories on isentropic surfaces to give, in some cases, greater 
error than from applying CAV trajectories to constant pressure sur- 
faces. This was not observed in the data collected in this study. 
This nay suggest that the assumption of no horizontal divergence is 
not too restrictive, and that possibly one of the other assumptions 
made by Rossby is more restrictive, e.g., assimiption of pxurely 
horizontal motion. 

If xve take a large constant value for horizontal velocity 
divergence, say 6 x 10~® per second, take a like vailue for the 
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shear term, assvme a time period of t^/elve hours, and use these 
values in equation (11), we would find that the value of at, 
say, 45° is reduced by 33, o or to the value of-f at latitude 33°. 
But, this is cin extreme case, and the effects of the divergence 
term are usually much less. However, since the assumption of 
does lead to error, it seems that some futiure re- 
search could be profitably spent introducing quantitative values 
for the divergence term and computation of a trajectory incorpor- 
ating this correction to the CAV trajectory for use in more 
accurate forecasting. 
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FREQUEN’CY 




ERROR 



Distribution of t^venty-four hour Aidnd direction 
error from (a) isentropic charts cind (b) from 
500 millibar charts 
Zrror: tens of degrees 

Figure 1. 
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FREQUENCY 




Distribution of forty-eight hotir wind direction 
error from (a) is en tropic charts and (b) from 
500 millibar charts 
Error: tens of degrees 

Figure 2. 
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Distribution of twenty-four hour ^riLnd speed 
error fron (a) isentropic charts and (b) 
from 500 millibar charts 
Error ; loiots 

Figure 3. 
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Distribution of forty-eight hour \d.nd speed 
error from (a) isentropic charts and (b) 
from 500 millibar charts 
Error; knots 



Figxire 4. 
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